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THERMODYNAMIC ANALYSIS OF AN AXIAL-FLOW 
LASER-POWERED THRUSTER 
by Albert F. Kascak 
Lewis Research Center 

SUMMARY 

A preliminary analysis of an axial -flow laser -powered thruster is performed. An 
axial -flow laser -powered thruster is a device in which laser energy is beamed axially 
through an exhaust nozzle into a cylindrical chamber. The laser energy is absorbed by 
small solid particles suspended in the propellant gas. These particles then heat the gas 
by conduction. The propellant flows axially from the chamber inlet toward the exhaust 
nozzle. 

The model analyzed assumes two inviscid flow regions; a region heated directly by 
the laser beam and a region between the chamber wall and the laser beam. The flow in 
each region is characterized by a bulk velocity and temperature, and there is a common 
pressxire for both regions. The propellant is assumed to be a perfect gas with negligible 
solid particle thermodynamic properties and a constant optical absorption cross section. 
The one-dimensional inviscid flow equations are solved for each region with the same 
initial conditions. 

The results show that, in general, the fluid in the laser -beam region heats up, ex- 
pands, and causes a pressure gradient which accelerates the less dense laser-beam 
fluid more than the more dense coaxial fluid. Fluid is convected from the laser-beam 
region to the coaxial region. As the laser-beam fluid is heated, it begins to radiate 
heat to the fluid in the coaxial region, which in turn becomes heated and expands. This 
heating results in a pressure gradient which causes fluid to be transferred from the 
slower moving coaxial region to the faster moving laser -beam region and thus causes 
the flxiid in the coaxial region to accelerate faster than the fluid in the laser -beam region. 

It is concluded that for a given inlet Mach number there is a maximum laser inten- 
sity which will choke the flow. This maximum intensity increases as the inlet Mach 
number decreases. If the laser beam is the same size as the thruster tube, the inlet 
Mach number must be less than 0.14 to achieve a factor of 10 increase in the propellant 
gas enthalpy. 

For a hydrogen gas seeded with 10 percent tungsten and entering the thruster tube 
at 10 atmospheres and 300 K and a laser beam that is one -half the area of the tube, a 
tube length between 20 and 50 centimeters is required to vaporize the tungsten particles. 



INTRODUCTION 


Several authors have considered the possibility of using high -power lasers as an 
energy source for rocket propulsion (e.g. , ref. 1). Energy is converted into a laser 
beam at some convenient location. This energy is transmitted (by the beam) to a remote 
location. At the remote location, the energy is absorbed in a thruster, where it is used 
to heat a propellant. The heated propellant is expanded through a nozzle and produces 
thrust. A conceptual system is illustrated in figure 1. System performance of laser- 
powered thrvister systems is described in references 1 and 2. 

An axial -flow thruster of the type suggested in reference 1 is shown in figure 2. 

The laser energy is beamed axially throt^h an exhaust nozzle into a cylindrical chamber . 
This energy is absorbed by small solid particles suspended in the propellant gas. These 
particles then heat the gas by conduction. The propellant flows axially from the cham- 
ber inlet toward the exhaust nozzle. Part of the propellant is heated by the laser beam 
and expands. The rest of the propellant does not "see" the laser beam and is heated by 
reradiation. This difference results in pressure gradients which cause flow profile re- 
distributions. The coupling of the laser energy to the propellant and the subsequent flow 
redistribution are the subjects of this report. One of the major questions to be an- 
swered is How large must the chamber be to absorb the laser energy? 

The approach taken is to assume two inviscid flow regions: the region heated by the 
laser beam and the region between the chamber wall and the laser beam. The flow is 
characterized by a bulk velocity and temperature in each region and a common pressure 
for both regions. The one -dimensional inviscid flow equations are solved for each re- 
gion with the same initial conditions. The heat transfer between the two regions is as- 
sumed to be by thermal radiation. 


SYMBOLS 

A area 

a absorption coefficient 

constant defined in eq, (8b) 

Cg constant defined in eq. (8c) 

Cg constant defined in eq. (8d) 

H convective energy 
h enthalpy 
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I laser -beam intensity 

M Mach number 

JK convective momentum 

m convective mass 

N constant defined in eq. (8e) 

P pressure 

Q heat generation rate 

T temperature 

u axial velocity 

z axial distance 

Of constant defined in eq. (12c) 

/3 constant defined in eq. (13c) 

y ratio of specific heat at constant pressure to specific heat at constant volume 
6 constant defined in eq. (15b) 

e constant defined in eq. (15c) 

rj constant defined in eq. (15d) 

^ ratio of coaxial flow area to laser -beam area 
p. density 

cr Steffan-Boltzmann constant 

T optical distance based on initial absorption coefficient 

Subscripts; 

B laser -beam region 

C coaxial flow region 

max quantity evaluated at condition of maximum laser -beam intensity 

0 inlet region 

Superscripts: 

~ quantity nondimensionalized by initial value of quantity 
’ derivative with respect to distance or optical distance 

* velocity or enthalpy associated with mass transferred into region for which equa 
tions are written 


ANALYSIS 


This analysis of an axial -flow laser-powered thruster assumed that the propellant 
is a mixture of solid particles and a perfect gas. The solid particles absorb the laser 
radiation and then conduct the heat to the surrounding gas. The heat-transfer coeffi- 
cient between the particles and the gas is assumed to be so large that there is negligible 
temperature difference between the particles and the gas. 

The model used in the analysis is shown in figure 2. The cylindrical chamber is 
divided into two regions. The first region (B) "sees” the laser beam. The second re- 
gion (C) is coaxial to the first and does not "see" the laser beam. The propellant flows 
axially from the chamber inlet toward the exhaust nozzle. At the inlet (0), the flow is 
imiform in both regions and has the same properties. 

The axial mass, momentum, and energy convection are defined, respectively, as 

(la) 

(lb) 

(l c) 

where p is the density, u is the axial velocity, A is the cross-sectional area, and h 
" is the enthalpy in either region B or C. For invLscid one -dimensional flow, the axial 
momentum and energy equation are written, respectively, as 

= -AP' + u*m' (2a) 

(2b) 

where P is the pressure, m’ is the mass transferred from one region into the region 
of interest, and Q is the heat generation rate, (The superscript * indicates the 
velocity or enthalpy that is transferred with the associated mass transfer, and the re- 
gion of interest is the region for which the equations are written.) 

The equation of state for a perfect gas can be written as 


H’ = QA + (h* + I u*nm' 


m = puA 
Ji = um 



P 




(3) 


where solid particle thermodynamic properties have been neglected, and T is temper- 
ature, The absorption coefficient for a gas seeded with small solid particles is propor- 
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tional to the mixture density: 


JL - P 
^0 PQ 

where it is assumed that the mass fraction of seed is small compared to that of the 
gas, and the seed particle radius remains constant with temperature and pressure. 

The heat source terms for the two regions are 

^ ^B " j 

Qc = ^^"B ^ (tb - Tc) 

where the radiation term between the two regions is based on the transparent approxima- 
tionj and all the heat radiated from region B is deposited in region C. (This formulation 
neglects conduction between the two regions.) The intensity of the laser beam Ig is 
given by the radiation absorption law: 

^B " ^B^B 

The inlet to the chamber is a permeable wall. Flow enters the wall from a reser- 
voir at a very low temperature (T ~ 0)'. As the flow passes throxigh the wall, it must 
pick up the laser energy fallii^ on the inside of the wall Agig q. Since the flow leaves 
the wall at a temperature Tq, the heat balance on the wall is 

^B%, 0 "" ^^B ^0 2 ^o) 

where Ag is a constant equal to the laser-beam cross-sectional area. 

The laser intensity at the exit of the thrust chamber is given. In order to make the 
problem an initial value problem, the laser intensity at the inlet must be known. The 
required inlet intensity can be calculated from equation (7) and an inlet wall temperature 
limit. The initial value problem can then be integrated to the exit intensity. The length 
for this integration is thus the length of the thrust chamber. 

The variables can be nondimensionalized by dividing the variable by its inlet value 
and denoting the nondimensional variable by a superscript bar. The axial distance is 
transformed to the nondimensional optical distance by 

T = a^z (8a) 

The prime will now refer to the derivative with respect to the optical distance. From the 
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( 4 ) 

(5a) 

(5b) 


perfect gas law and thermodynamics 


^0 _ y - 1 = 

Pq\ y 

where y is the ratio of specific heats. The inlet Mach number Mq is 

1 ^0 , Cg 

0" 1 K 

r - 1 hp 1 (y . 1) 

2 

At the inlet the ratio of the radiated to convected heat transfer is 

C _ °^0 (8d) 

where N is the radiative to convective number, defined as 


aT’ 

N= 1- (8e) 


(8b) 


(8c) 


Using the constants defined in equation (8) and the previously defined nondimension- 
al variables gives for the momentum and energy , equations 


and 


Jl' 


2C. 


— AP' + u* m’ 


(9 a) 


H’ = 


QA 


(1 + CglaQPQUphQ 


h* + 

1 + C, 


(9b) 


where 


Q 


B 




= a 






(9c) 
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Q. 




= 4aC,^ 


(9d) 


(Eq. (9c) used the initial condition, eq. (7).) With the nondimensional variables defined, 
the equation of state (eq. (3)) becomes 


-= E = £ 

^ f h 


(10a) 


and from equation (3) the nondimensional absorption coefficient becomes equal to the 
nondimensional density ; 


a = p 

The radiation absorption law (eq. (6)) then becomes 


(10b) 


T ' P 

B _:r- _ B 

= 


(11) 


Using the definition of the nondimensional axial momentum convection resvilts in 

ui' = um (12a) 

and equation (9a), the derivative of the nondimensional axial velocity, becomes 

(12b) 


u 


C — -L. 

'^l P' n? 

+ a — 


m 


^^2 ^ 


where a is defined as 


a - 


u* - u 


(12c) 


u 


Using the definition of the nondimensional axial energy convection gives 


1 + C, 


(13a) 


and equations (9b) and (12b) ; the derivative of the nondimensional enthalpy becomes 

— 2 — + 
h puh aQpQUQhQ ph m 


(13b) 
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where p is defined as 




Ti* - h + C2(u* - u)' 


(13c) 


Using the definition of the nondimensional axial mass convection gives 

in = puA (14a) 

and the equation of state (10a); the derivative of the nondimensional axial mass convec- 
tion becomes 


nT + IL -Jll + — 

mu P h A . 

Substituting equations (12b) and (13b) into equation (14b) and rearranging yield 

« in' ^ P’ . „ 
m P 

where 


6 = 1- O' + /3 


e = 1 - Cj - 





^ Q 

puh agPo^O^O 


(14b) 

(15a) 

(15b) 

(15c) 

(15d) 


The mass leaving one region must be equal to the mass entering the other; 
therefore, 

nig = -mQ (16) 

It is assumed tiiat the pressures in both regions are equal, that is, there are no radial 
pressure gradients; 

P = Pb = Pc (17) 

Substituting equations (10a), (14a), and (15a) into equations (12b), (13b), and (16) yields 
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a€ 

.6 


(18a) 


Hi 

u 


2C, 



P* , ay 
P 6 


=(fii +0,) El t a +1- 

h \6 VP 6 A 

P' ^ ^B^B 
^ ^B^B ^ ^ "C^C 

^B®B ^C^C 


(18b) 


(18c) 


where equations (18a) and (18b) are applicable in both region B and region C. Substi 
tuting eqxiations (10a) and (14a) into (15a) and evaluating it in region B yield 



^B®B 



(19) 


If < 0, mass is being convected from region B to region C, and the associated ve- 
locity and enthalpy are 


U* =Ug 

(20a) 

h* =hg 

(20b) 


If > 0, mass is being convected from region C to region B, and the associated ve- 
locity and enthalpy are 


II 

(21a) 

h* =h(. 

(21b) 


Equations (11) and (18), when applied to both regions, form a complete set of six 
differential equations in six unknowns: Ug, u^, Hg, P, and I^. These equations 
can be solved by using a predictor -corrector nvunerical technique. Initially the vari- 
ables are normalized to 1. Auxiliary equations used are (9) , (10), (12), (13), (15), (19), 
and either (20) or (21). The gamma of the gas, the initial Mach number, the initial 
radiative to convective number, and the ratio of the wall area to the beam must be given. 
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A special case arises when the wall area is equal to the cross-sectional area of 
the beam. For this case there is no flow outside the beam, and the wall temperature is 
equal to the beam gas temperature: 

= 0 (22a) 

Ag = 1 (22b) 

Tg = T^ (22c) 


The mass convection in the beam is a constant: 


Therefore, 


raB = l 




(23a) 


(23b) 


The momentum equation (9a) can be integrated and when solved for the pressure becomes 

2C, 


P = <"B-« 


(24a) 


Using equations (23b) and (24a) gives the enthalpy: 

2C 


^B “ “b 


(Ub - 1) 


The heat source term (9c) becomes 


(24b) 


Q 


B 




- = (1 + C2)alg - (1 + C,)!^ 


2':"B 


Substituting this into the energy equation and integrating yield 

hg + C2^| - (1 + C 2 ) = (1 + C2)(rB - 1) 


(25a) 


(25b) 


Substituting equation (24b) into (25b) and solving for the intensity give 
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DISCUSSION 

A model of a coaxial flow laser -powered thruster is shown in figvire 2. The cylin- 
drical chamber is divided into two coaxial regions. The first (inner) region is heated 
directly by the laser beam. The propellant flows axially from the chamber inlet toward 
the exhaust nozzle. At the inlet the flow is uniform in the two regions, both having the 
same temperature, pressure, and velocity. The fluid in the center is heated and ex- 
pands and thus causes pressure gradients. These pressure gradients accelerate both the 
heated and the nonheated fluid. As the flow accelerates, energy and mass are trans- 
ferred between the two regions. The laser is assumed to be beamed through the nozzle. 
It is further assumed that the beam intensity in the propellant decreases exponentially 
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with the distance from the nozzle. The region of interest is shown in figure 2. 

The first case studied is one in which tiie laser beam occupies the entire cylindri- 
cal volume. There is no coaxial flow. This condition is equivalent to heating a tube 
with nonfrictional flow. It has been analyzed in many texts on advanced fluid flow, such 
as reference 3. These analyses showed that there is a maximum amount of heat that 
can be deposited in the gas in the tube. This condition is expressed in equation (28), in 
which the laser-beam intensity is maximized. With this value of the axial velocity 
known, the corresponding variables are obtained from equations (24), (26), and (27) and 
are plotted in figure 3 . 

Figure 3 shows nondimensional axial velocity, enthalpy, laser intensity, pressure, 
and optical distance plotted against the initial Mach number for the condition of maxi- 
mum laser -beam intensity. For this case the laser beam fills the thruster tube, and 
there is no coaxial flow. Values for a monatomic gas, such as dissociated hydrogen 
(y = 1.667), are shown in figure 3(a); and values for a diatomic gas, such as undisso- 
ciated hydrogen (y = 1.4), are shown in figure 3(b). 

Figure 3 shows that to increase the maximum laser -beam intensity, and thus the 
amount of energy deposited in the propellant gas, the initial Mach number must be de- 
creased. If the propellant gas is heated to at least 10 times the initial enthalpy, the 
initial Mach number must be less than 0.14, and the minimum optical distance must be 
greater than 13. 

An example of the detailed axial profiles which result in this maximum laser -beam 
intensity is shown in figure 4. This figure shows the Mach number and nondimensional 
axial velocity, enthalpy, laser intensity, and pressure plotted against optical distance 
(the initial absorption coefficient times the distance from the inlet). The initial Mach 
number is 0.01, and the laser fills the thruster tube (there is no coaxial flow). Figure 
4(a) is for a monatomic gas, and figure 4(b) is for a diatomic gas. Over a large portion 
of the range of optical distance there is very little difference between a monatomic and a 
diatomic gas; this is the region where the Mach number is low. As the Mach number 
approaches 1, the flow begins to choke. The slopes of all the curves in figure 4 ap- 
proach plus or minus infinity. As the Mach number increases above 1, the laser-beam 
intensity decreases. Since the laser-beam intensity must increase, this portion of the 
curve is not physically possible. 

In any real system, the laser -beam cross-sectional area would be less than the 
cross-sectional area of the thruster tube. (Besides the obvioxis problems of misaline- 
ment and heat transfer to the tube wall, the previous discussion indicated a need for a 
convergent -divergent nozzle to make the flow supersonic. This nozzle shadows a por- 
tion of the cross-sectional area of the tube.) Figure 5 presents results for a case in 
which the cross-sectional area of the laser beam is one-half the cross-sectional area 
of the thruster tube. In this figure the calculated Mach number and nondimensional 
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axial velocity, enthalpy, mass convection, and intensity in the laser-beam region; Mach 
number and nondimensional axial velocity, enthalpy, and mass convection in the coaxial 
region; and nondimensional pressure are plotted against optical distance from the inlet. 
At the inlet the initial pressure was 10 atmospheres, and the initial temperature was 
300 K. Figures 5(a) and (c) represent undissociated hydrogen, and figures 5(b) and (d) 
represent dissociated hydrogen. Figures 5(a) and (b) present data for an initial Mach 
number of 0. 01, and figures 5(c) and (d) for an initial Mach number of 0.005. 

In general, the fluid in the laser -beam region heats up and expands. This expan- 
sion causes a pressm*e gradient which accelerates the fluid in the laser -beam region 
i^d to a lesser extent the fluid in the coaxial region, (The latter fluid is accelerated 
less because it is more dense.) Some of the fluid in the laser-beam region is convected 
into the coaxial region, and thus the enthalpy in the coaxial region is raised. As the 
enthalpy in the laser -beam region is raised, the region begins to radiate energy to the 
coaxial region. The coaxial region then heats up and expands. This expansion causes 
a pressure gradient which accelerates the fluid in the coaxial region and to a lesser ex- 
tent the fluid in the laser region. (Some of the more slowly moving fluid in the coaxial 
region is convected back into the laser region and reduces acceleration in the laser 
region.) 

If the solid seed particles in the propellant gas were tungsten, the particles would 
vaporize between 5000 and 6000 K, This behavior corresponds to an enthalpy increase 
of a factor of about 20. An enthalpy increase of a factor of 20 in the laser -beam region 
corresponds to an optical distance of about 14 for an initial Mach number of 0.01 and an 
optical distance of about 21 for an initial Mach number of 0.005. For a seed mass frac- 
tion of about 10 percent and a seed absorption cross section 10 000 square centimeters 
per gram, the corresponding thruster tube length is between 20 and 50 centimeters. 


CONCLUSIONS 

The followii^ conclusions were drawn from a thermodynamic analysis of a two- 
region model of an axial-flow laser-powered thruster: 

1. There is a maximum laser intensity which will choke the flow in the thruster 

tube. 

2. If the inlet Mach number is reduced, the maximum laser intensity can be in- 
creased. 

3. If the laser beam is the same size as the thruster tube, the inlet Mach number 
must be less than 0. 14 to achieve a factor of 10 increase in the propellant gas enthalpy. 

4. For hydrogen gas seeded with 10 percent tungsten and enterii^ the thruster 


13 



tube at 10 atmospheres and 300 K and a laser beam that is one -half the area of the tube, 
a tube lei^th between 20 and 50 centimeters is required to vaporize the tungsten parti- 
cles. 

Lewis Research Center , 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 1, 1974, 

502-04. 
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Figure 1. - Laser rocket system. 


Inlet region, 0 

Uniform propellant flow 
seeded with small solid 
particles 



Figure 2. - Model of coaxial flow laser-powered thruster. 


Nondimensional axial velocity, enthalpy, laser intensity, pressure, or optical distance 




Figure 3. - Nondimensional axial velocity, enthalpy, laser intensity, pres- 
sure, and optical distance as functions of initial Mach number for condi- 
tion of maximum laser intensity. Laser fills thruster tube (no coaxial 
flow). 




Mach number or nondimensional axial velocity, enthalpy, laser intensity, or pressure 




Figure 4. - Mach number and nondimensional axial velocity, enthalpy, laser 
intensity, and pressure as functions' of optical distance for initial Mach 
number of 0.01. Laser fills thruster tube (no coaxial flowl. 



Laser-beam region Mach number or nondimensional axial velocity, enthalpy, mass convection, or intensity; coaxial region Mach number or nondimen- 

sional axial velocity, enthalpy, or mass convection; or nondimensional pressure 



(a) Specific-heat ratio, L667; molecular weight, 1; initial Mach number. 



(b) Specific-heat ratio, 1.4; molecular weight, 2; initial Mach number, 

a 01. 


Figure 5. - Laser-beam region Mach number and nondimensional axial 
velocity, enthalpy, mass convection, and intensity; coaxial region 
Mach number and nondimensional axial velocity, enthalpy, and mass 
convection; and nondimensional pressure as functions of optical dis- 
tance. Laser beam fills one-half of thruster tube; initial pressure, 

10 atmospheres; initial temperature, 300 K. 



Laser-beam region Mach number or nondimensional axial velocity, enthalpy, mass convection, or intensity; coaxial region Mach number or nondimen- 

sional axial velocity, enthalpy, or mass convection; or nondimensional pressure 



(c) Specific-heat ratio, 1.667; molecular weight, 1; initial Mach number, 
a 005. 



Optical distance, t 


(d) Specific-heat ratio, 1.4; molecular weight, 2; initial Mach number, 
a 005. \ 

Figures. - Concluded. 


NASA-Langley, 1974 


E-7891 


19 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



POSTMASTER ; 


It Undeliverable (Section 158 
Postal Manual) Do Not Return 


"T/>e aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemiruttion 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 





